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Abstract 

We study the m = 3 bootstrap percolation model on a cubic lattice, 
using Monte Carlo simulation and finite-size scaling techniques. In boot- 
strap percolation, sites on a lattice are considered occupied (present) or 
vacant (absent) with probability p or 1 — p, respectively. Occupied sites 
with less than m occupied first-neighbours are then rendered unoccupied; 
this culling process is repeated until a stable configuration is reached. We 
evaluate the percolation critical probability, p c , and both scaling powers, 
y v and yh, and, contrarily to previous calculations, our results indicate 
that the model belongs to the same universality class as usual percolation 
(i.e., m = 0). The critical spanning probability, R(p c ), is also numeri- 
cally studied, for systems with linear sizes ranging from L = 32 up to 
L = 480: the value we found, R(p c ) ~ 0.270 ± 0.005, is the same as for 
usual percolation with free boundary conditions. 
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1 Introduction 

In many physical phenomena percolation effects play an important role Q. 
Particularly, some dilute magnets are well described, in what concerns magnetic 
phase transitions, by uncorrelated diluted models ||. In these models, magnetic 
sites on a lattice are randomly replaced by non-magnetic ones, and a bond links 
each pair of occupied (magnetic) first neighbours. At zero temperature, the 
problem is purely geometrical and is described in the following way. Sites on a 
lattice are randomly occupied with probability p, while all bonds are considered 
to be present. A cluster is then defined as a collection of present sites which are 
connected to each other through steps between occupied first-neighbours. As p 
increases, an infinite cluster appears for the first time at a critical probability p c , 
which is lattice dependent. In analogy with thermal critical phenomena, some 
quantities are singular at the critical point, following a power-law behavior near 
p c Q. Typical examples are the probability P(p) that a site belongs to the 
infinite cluster, which behaves as P(p) ~ (p — p c ) near p c , and the correlation 
length £(p), which diverges at p c , according to £(p) ~ \p — p c \~ v ■ 
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Nevertheless, some systems were found to be better described by correlated 
percolation models, where the presence of sites (or bonds) depends also on 
their neighbourhood. Typical examples of correlated percolation models are 
bootstrap percolation B and site-bond correlated percolation Q] . The physical 
motivation for the introduction of the former model comes from diluted magnetic 
systems where competition between exchange (which favour a magnetic ground 
state) and crystal-field (which leads to a non-magnetic ground state) interactions 
takes place. To mimic this competition at zero temperature, the bootstrap 
percolation model was introduced || [|: in this model, sites on a lattice are 
randomly occupied with probability p but only those with at least m occupied 
first-neighbours remain occupied. In the stable final configuration, all occupied 
sites have at least m occupied first-neighbours or the whole lattice is empty. 
As our purpose in this paper is to study one implementation of the bootstrap 
percolation model, we briefly review some of its properties in what follows (for 
a thourough discussion of the results up to 1990, see Reference ||). 

The m — case regains the usual (uncorrelated) site percolation model, 
where the transition is continuous and p c < 1 in two and three dimensions. 
More specifically, the most precise evaluations of the critical exponents for three- 
dimensional lattices are v = l/y p = 0.875 ± 0.008 and (3 = 0.412 ± 0.010 from 
simulation § and v = l/y p = 0.872 ± 0.070 and f3 = 0.405 ± 0.025 from 
series In the m — 1 case, only isolated sites are removed by the culling 
process: these sites do not contribute to the critical behaviour and p c and all 
critical exponents remain the same as for usual percolation. For m = 2, on 
the other hand, isolated clusters with two sites are eliminated, as well as some 
dangling structures of more compact clusters. The elimination of these dangling 
structures, however, does not break the infinite cluster (whenever present) and, 
therefore, the critical probability is the same as for usual percolation. Moreover, 
as the exponent v is also connected to the formation of this infinite cluster, its 
value is the same for m = 2 and m = 0. In what concerns "field" exponents 
(/3, for instance), previous results for two-dimensional systems indicate a higher 
value of this exponent for m — 2 than for usual percolation Nevertheless, it 
has been shown later, through simulations on bigger lattices in two dimensions 
[ [To| and general arguments applied to both two and three dimensions ]ll) , that 
the exponent [3 is the same for uncorrelated and m = 2 bootstrap percolation 
models. 

Let us now turn our atention to higher values of m. It is generally believed 
that, for any value of m where only infinite clusters can survive the culling 
process, p c = 1. This is indeed the case for m = 2d — 1 on hypercubic lattices 
(d is the dimension of the lattice) |l^] , for m = 4 on cubic |l^] and triangular 
[|| lattices and for m = 5 on the triangular lattice jljS. Moreover, it has been 
shown that, for these cases, the usual finite-size scaling relation does not hold. 
This finite-size scaling predicts that, if a suitable definition of a finite-lattice 
"critical" point, p av , is made, this point will approach the critical point in the 
thermodynamic limit, p c , as: 

p c -Pav~ L- 1 '", (1) 
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where L is the linear size of the finite system, such that L ^> 1, and v is the 
usual critical exponent. This finite-size behaviour is indeed observed for m < 2, 
but fails for high values of m. For m — 2d — 1 on hypercubic lattices, it has 
been proven that p c —p a v ~ l/(l°g d 1 L) @, e.g., proportional to l/log(logL) 
for d = 3 dimensions. Also for m = 4 on the cubic lattice, the correct finite-size 
behaviour is p c — p av ~ 1/ log (log L) , with p c =\ |0 . These results for high m 
have been conjectured or tested in numerical simulations Jl6| . 

The results stated in the previous paragraphs do not apply to the m = 3 
case on the cubic lattice. For this model, one expects p c to be above the value 
for uncorrected percolation (where p c = 0.311605 ± 0.000005 @), since the 
infinite cluster for usual percolation at p c is not stable with respect to the culling 
process for m = 3. On the other hand, since finite clusters are still stable for this 
value of m on the cubic lattice, it is expected that p c < 1. Numerical simulations 
confirmed this scenario |]l3[ [ll|, although the relative small sizes used in those 
works indicate that the values might not be precise (we will return to this point 
later). In what concerns the critical exponents, previous results indicate that 
the exponent v is the same as for usual percolation JOI but is higher than 
its uncorrelated counterpart jl3|, However, in neither of these works a 

extrapolation to the thermodynamic limit is attempted, leaving the possibility 
that finite-size effects might be the reason for the discrepancy in the values of (3. 
This possibility was first proven to be right, in the context of two-dimensional 
bootstrap percolation models, in References jl0| and later confirmed for m = 3 
on the triangular lattice pl|. The possibility of a new universality class for 
bootsrap percolation with m = 3 on the cubic lattice is the problem we address 
in this work. We resort to numerical simulation methods, which, together with 
finite-size scaling analysis, allowed us to obtain more precise values for the 
critical parameters. We also study the so-called critical spanning probability, 
i.e., the probability that a given lattice has a cluster connecting its boundaries at 
criticality |2(], |2l| . This quantity shows some degree of universality: it depends 
on the dimension and shape of the system and on the specific boundary condition 
but not on the lattice type (simple cubic of f.c.c, for instance) and on the 
particular kind of percolation (site or bond) 0, |2l], 22 . It is then interesting 
to see whether it remains invariant for percolation models such that long-range 
correlation is involved, like bootstrap percolation. 

The remainder of the paper is organized as follows. In the next section we 
present the method and discuss some technical details, as well as the results for 
the critical parameters. In Section^ we discuss some previous results concerning 
the critical spanning probability for percolation and our results for bootstrap 
percolation. Finally, we summarize our results in the last section. 



2 Method and Results 

The method we use is connected to real-space renormalization-group and finite- 
size scaling procedures j^3|. The approach needs that precise values of the 
physical quantities are available; these are only obtained for high values of the 
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linear system size L. We then studied finite systems of size L 3 , with 32 < L < 
480; from the results for L ^> 1, it is possible to use finite-size scaling techniques 
to extrapolate to the thermodynamic limit (L = oo). 

The critical probability p c and the critical exponent v are calculated as 
follows. For a lattice of size L, we occupy each site with probability p, apply 
the bootstrap condition and test the lattice for percolation (here we define 
percolation as the presence of a cluster which connects the bottom and top 
planes of the finite cubic lattice; we discuss this and other technical points 
below). Our finite-size estimate of the critical probability p* , is taken as the 
value of p at which the cell percolates for the first time, when p is increased 
from zero; this procedure is made for TV different runs (which correspond to TV 
different seeds to the random number generator). Each run leads to a different 
value of p* , since the lattice is finite. We take the average of the TV values 
of p* as our estimate of p av (see Section [l]). It is then assumed that p av will 
approach p c as given by Equation [j]. Moreover, it is possible to calculate the 
width a = y/ < p* 2 > -~Pa V , which behaves as (< p* > stands for the average 
of p* 2 over the TV realizations) : 

a-L- 1 ^. (2) 

From the previous equation and a log — log plot of a X L, it is then possible 
to obtain the value of the critical exponent v. Tha data is depicted in Figure 
|l|; from the slope of the straight line we obtain v = 0.89 ± 0.04. We compare 
this value with other evaluations in Table [l]; within the numerical precision, 
this exponent is the same for bootstrap percolation with m = 3 and ordinary 
percolation on the cubic lattice and agrees with previous evaluations of v for 
bootstrap percolation with m = 3 on the cubic lattice. Note from the graph 
that the straight line regime is achieved for L > 128, while for uncorrelated 
percolation this regime sets in for smaller values of L. This is expected, since 
finite-size effects are stronger for correlated percolation problems than for their 
uncorrelated counterparts || [24J. Therefore, we neglect the data for L < 96 
and used only lattices with 128 < L < 480 in our linear regression analysis. 

In order to calculate the critical threshold, we resort to Equations [l] and ||, 
which imply that, for L ^> 1: 

Pc - Pav ~ er. (3) 

This is a convenient way to calculate p c , since it does not depend on the value of 
u; such dependence would appear if Equation |l| was used. It is shown in Figure 
H a plot of Pav x Pc is given by the linear coeficient and the value obtained 
is p c — 0.57256 ± 0.00006. This value is slightly above the one calculated in 
Reference [Q, in which the value of L varied from 10 to 110. If we use the 
same range in our calculation the extrapolated value of p c is consistent with the 
result of g§ (see Table |l|). 

Let us know discuss some technical points. An important quantity is the 
value of p at which the lattice percolates for the first time, when p is increased 
from zero. One needs to define what "percolate" means for a finite lattice. In 
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Figure 1: Log-log plot a x L, for 32 < L < 480. The linear size L is always a 
multiple of 32. The dotted line is a linear fit of the data for 128 < L < 480 only 
(see text). 
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Figure 2: Graph of p av x a, for 32 < L < 480. The linear size L is always a 
multiple of 32. The dotted line is a linear fit of the data for 128 < L < 480 
only (see text). Although not depicted in the figure, there are error bars in the 
values of a, which were taken into account in the evaluation of p c - 
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Table 1: Critical parameters for m — 3 bootstrap percolation and comparison 
with previous results for this model and for usual percolation. 





m=3 bootstrap percolation 


usual percolation 


Pc 


0.57256±0.00006 a 






0.5717±0.0005 fc 


0.311605±0.000005 d 




0.568±0.002 c 
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0.89±0.04° 


0.875±0.008 e 




0.876±0.010 b 


0.872±0.070 / 





0.37±0.03 a 


0.41±0.02 e 




0.6±0.1 fc 


oAo±om f 




0.82±0.04 c 
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this work, we used the rule called TZi in Reference |23j], i.e., a lattice percolates 
if, after the bootstrap culling process, there is a path of present sites which links 
the boundaries of the lattice in a fixed direction (vertical, say). There are other 
possible definitions (see |23| ) and it is expected that all of them lead to the same 
value of p c in the thermodynamic limit. Note, however, that the value of the 
critical spanning probability does depend on this definition, as we will discuss 
im Section |^. The numerical procedure we used to test for percolation is the 
Hoshen-Kopelman algorithm pq ): for usual percolation, it requires the storage 
of only one plane. For bootstrap percolation, on the other hand, the bootstrap 
iteration needs the storage of the whole lattice, due to correlation effects. To 
cope with this drawback, we store the lattice in bits, instead of words: this 
saves memory and time, since the updates connected to the bootstrap rule can 
be made in parallel for a set of 32 sites |2(| . To define all six first-neighbours of 
sites at the boundaries of the lattice, periodic boundary conditions are used. It 
is expected that the boundary condition does not affect the values of the critical 
parameters in the thermodynamic limit, since it is a "surface" effect. Finally, 
let us mention that the number of realizations M (see Section |l|) varied from 
12000 for L = 32 to 640 for L = 480; the errors were calculated as three times 
the standard deviation for subsets of the total number of realizations. 

To have access to the "magnetic" scaling power y h , one possible way is to 
define a "ghost"- site, which is linked to all sites of the lattice with probability 
h pjj| . Within a real-space renormalization group framework, it is possible to 
calculate the eigenvalue through \ h =< n > /p c , where < n > is the average 
number of occupied sites linked to one of the boundary planes of the finite cubic 
lattice at the critical threshold, p = p c . The value of p c is obtained as explained 
above and is then used to calculate < n > and hence \h(L), averaging over M 
configurations. The value of Ai varied from ~ 320000, for the smaller lattices, 
to ~ 20000, for the bigger ones. We use two procedures to obtain yh in the 
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thermodynamic limit. The first one is based on the fact that, for L — > oo: 

\ h {L) = D>\ (4) 

with yh independent of L; this equation leads to a straight line in a log — log 
plot of Xh x L, for L>1. As depicted in Figure g| this is indeed the case for 
128 < L < 480. From a linear fitting of the data for this range of L, and resorting 
to the scaling relation j3 = (d — yh)v, we obtained the value (3 — 0.37 ±0.03 (see 
Table |l|). Alternatively, we could take into account correction- to-scaling terms, 
through: 

X h = I?* (1 + B/L) , (5) 

From this equation, we see that local slopes of log x logL provide estimates 
of yh{L), which, when extrapolated to L — » oo, leads to an evaluation of yh in 
the thermodynamic limit. We have applied this procedure for 32 < L < 480, 
using three consecutive points to caudate the local slopes and extrapolating to 
L — > oo through a yh{L) xl/L graph. The value thus obtained for /3 is the same 
as for the first procedure. In Table |l| we see that our estimate of j3 disagrees with 
the values calculated in References [[l3] and In evaluating this exponent, 
the former uses lattices of linear size L — 35, while the latter uses L = 80; 
neither attempted to make an extrapolation to the thermodynamic limit. From 
Table |, we can infer that our estimate of (3 for the bootstrap model we study 
is, within the error bars, equal to the corresponding value of this exponent 
for usual (uncorrelated) percolation. Since we have already seen that v is also 
the same for both models, we can draw the conclusion that usual percolation 
and bootstrap percolation with m = 3 on the cubic lattice belong to the same 



universality class. This result contradicts References |13 and [p|; we believe 
this is caused by the small lattices used in those works. 



3 Critical spanning probability 

It has been established some time ago that the critical spanning probability, 
R{p c ), defined as the probability of spanning a lattice at the critical point, 
shows some degree of universality |^0|, ||]J . More precisely, this quantity does not 
depend on the lattice type and on the kind of percolation. This result contradicts 
the assumption made on early applications of real-space renormalization group 
procedures to percolation. In those, it was assumed that the critical spanning 
probability is equal to the critical percolation threshold, p c and, therefore, lattice 
dependent |f23f . However, later numerical tests confirmed and expanded the 
universality proposal [|l7[ pl[ p2[ . 

Nevertheles, no study has been made on correlated models, to the best of 
our knowledge. While it is expected that short-range correlations do not change 
the universality scenario it is not clear whether R(p c ) changes when long- 
range correlations are introduced. A convenient model to test these possibilities 
is the bootstrap percolation one. While the behaviour of R{p c ) is trivial for 
the cases where p c — 1, for m — 3 on the simple cubic lattice the presence of 
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Figure 3: Graph of log(X h (L)) x log(i), for 32 < L < 480. The linear size 
L is always a multiple of 32. The dotted line is a linear fit of the data for 
128 < L < 480 only (see text) and yh is the angular coefficient of this line. The 
error bars are smaller than the data points. 



correlation may lead to a non trivial behavior. We study this possibility using 
numerical simulation on simple cubic lattices of size L 3 , with 32 < L < 480. The 
programs and algorithms used are essentially the same as the ones described in 
the previous section. The basic procedure is to generate M. independent runs 
for each lattice size and compute the fraction of those which percolate after the 
bootstrap condition is used and a stable configuration is reached. The values of 
M. are the same as those used in the calculation of yh (see previous section). 

The results are depicted in Figure^: we can infer that R(p c ) — 0.270±0.005, 
where the error bar is a rough estimate. There are two previous calculation 
of R(p c ) for the uncorrelated site percolation on the cubic lattice with free 
boundary conditions: they lead to R(p c ) = 0.265 ± 0.005 ||] and R{p c ) = 0.28 
| jl7| . Our value agrees, within the numerical error, with the first one but the 



result of Reference 17 cannot be ruled out. It is then reasonable to infer 
that usual percolation and to ~ 3 bootstrap percolation on the cubic lattice 
belong to the same universality class, also in what regards the critical spanning 
probability. 



4 Summary 

We calculate, using numerical simulations and finite-size scaling techniques, the 
critical parameters of the bootstrap percolation model with to = 3 on the cubic 
lattice, using finite lattices of size L 3 , with 32 < L < 480. Our evaluations of v 
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Figure 4: Graph of R(p c ) x log(L), for 32 < L < 480. The linear size L is always 
a multiple of 32. The dotted line is a guide to the eye and indicates the value 
0.270, taken as R{p c ) in the limit L — > oo. 



and (3 strongly support the conclusion that usual percolation and the model we 
study belong to the same universality class. This result disagrees with previous 
calculations Jl3, 18 ; we believe this is due to the small sizes used in those works. 
To support this assumption, we note that the finite-size scaling assumptions hold 
for lattices of linear size L > 128, which are higher then the sizes studied in 
previous works. 

The critical spanning probability, R(p c ), is also calculated. It has been 
shown that this quantity shows some degree of universality, but, to the best of 
our knowledge, no study concerning correlated percolation models have been 
done so far. Our result for bootstrap percolation with m = 3 on the cubic 
lattice, R{pc) — 0.270 ±0.005, is, within the numerical accuracy, the same value 
as for usual percolation with free boundary conditions. Therefore, we can infer 
that R{p c ) is not sensitive to short range correlations and even to some long 
range correlations, like the one studied in this paper. 
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of this work and for a critical reading of the manuscript. 
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